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Abstract 

The conventional paradigm of charge order for Lai-xCa^MnOs for x = 0.5 has been challenged 
recently by a Zener polaron picture emerging from experiments and theoretical calculations. The 
effective low energy Hamiltonian for the magnetic degrees of freedom has been found to be a 
cubic Heisenberg model, with ferromagnetic nearest neighbor and frustrating antiferromagnetic 
next nearest neighbor interactions in the planes, and antiferromagnetic interaction between planes. 
With linear spin wave theory and diagonalization of small clusters up to 27 sites we find that 
the behavior of the model interpolates between the A and CE-type magnetic structures when a 
frustrating intraplanar interaction is tuned. The values of the interactions calculated by ab initio 
methods indicate a possible non-bipartite picture of polaron ordering differing from the conventional 
one. 

PACS numbers: 75.10.Jm Quantized Spin Models, 75.30.Ds Spin waves, 75.47.Lx Manganites, 75.50.Ee 
Antiferromagnetics 
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I. INTRODUCTION 



Charge, orbital and spin order in doped manganites (Ai_ x B x Mn03, where A is a trivalent 
ion and B a divalent ion) with x > 0.4 is generally discussed in terms of the Goodenough 
model ]J and the double exchange Hamiltonian [2]. In that picture, charge ordering (CO) 
and orbital ordering (00) is associated with 3d populations of Mn ions, which may be d 3 or 
d 4 . In Lao.sCao.sMnOs, for example, it is believed that there are equal numbers of Mn ions 
with d 3 and d 4 populations in the CO state |3j . However, recent experimental and theoretical 
investigations point to a more complex picture, and there is evidence that a Zener polaron 
or similar model, which is quite distinct from the original Goodenough model, describes the 
CO state of Lao.sCao.sMnOs and Pr o .6oCa o .4oMnO 3 0, 0, E|. Here a frustrated Heisenberg 
Hamiltonian is used to model the Zener polaron state of half-doped manganites using spin 
wave theory and exact diagonalization of clusters with periodic boundary conditions. The 
parameter space for the Hamiltonian is derived from ab initio calculations on the Zener 
polaron state of Lao.sCao.sMnOs 0. In particular, we report energetics, spin waves, mag- 
netization, spin correlation functions, the spin contribution to heat capacity and magnetic 
susceptibility for this model. 

In the Zener polaron picture pairs of Mn ions are tightly bound into ferromagnetic (FM) 
dimers which interact relatively weakly with each other |4j,|5|. The Mn ion 3d population is d 4 
on all Mn ions and electrons are transferred from O 2- ions located between Mn ions in dimers, 
to form Mn-0~-Mn polarons. Electron transfer is necessary to satisfy electron counting. 
The dimers (Zener polarons) have spin 7/2. Their electronic structure is reflected in the 
Pm crystal structure of Pr .6oCa .4oMn0 3 determined by neutron scattering from a single 
crystal [4j. The P2i/m structure for Lao.sCao.sMnOs determined by neutron scattering]^ 
from a powder is different in several respects. The crystal structure of Lao.sCao.sMnOs has 
been calculated by using ab initio methods to minimize the crystal total energy [7| and the 
resulting structure resembles that for Pr .6oCao.4oMnOs The Zener polaron picture for 
half-doped manganites may therefo re apply more generally In the frustrated Heisenberg 
Hamiltonian used in this work, polarons are treated as single magnetic units with spin 7/2. 
This is appropriate for Mn ions tightly bound into FM dimers when the internal FM coupling 
greatly exceeds the inter-polaron coupling. Ab initio calculations on Lao.sCao.sMnOs show 
that magnetic coupling within polarons is strong and ferromagnetic (J ~ 200 meV), while 
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coupling between polarons along their zig-zag chains is also FM but much weaker (J ~ 10 
meV) p. There are both antiferromagnetic (AF) and FM couplings of the same order of 
magnitude between polarons in neighboring chains so that interchain coupling is strongly 
frustrated. 

Zener polarons and exchange couplings in Lao.sCao.sMnOa are shown schematically in Fig. 
together with the magnetic unit cell. Intrachain FM coupling is labelled FM1. Competing 
AF and FM intraplanar, interchain couplings are labelled AF1 and FM2; there is also an 
AF interplanar coupling, AF2, which is not shown. Ab initio calculations on LaMnOs show 
that exchange coupling constants in that compound depend on 00 and that not only the 
magnitudes, but also the signs of exchange constants can change when 00 changes |9|- In 
both LaMnOa and Lao.sCao.sMnOa, AF couplings are found between neighboring Mn ions 
when their filled e g orbitals are both oriented perpendicular to a Mn-O-Mn axis, while FM 
couplings are found between Mn ions when their filled e g orbitals are oriented perpendicular 
to each other, one being oriented along the Mn-O-Mn axis and the other perpendicular to 
it. FM exchange constants FM1 and FM2 in Lao.sCao.sMnOs were found to be similar in 
magnitude by ab initio calculation (—12 and —14 meV, respectively) they are assumed 
to have the same value, Ji, in this work. AF1 and AF2 exchange constants were found to 
be 5 and 8 meV, respectively, and correspond to J2 and J3 in the model described below. It 
is noted that Unrestricted Hartree-Fock (UHF) calculations used to obtain these exchange 
constants underestimate the magnitude of the AF exchange constants in LaMnOs, while 
the FM exchange constant is in agreement with the experimental value [9]. Inspection of 
Fig. E^a) shows that each polaron is ferromagnetically coupled to four neighbors, two on 
each side, and antiferromagnetically coupled to two neighbors in the same plane (Fig. Cfb)), 
along the diagonal direction; this arrangement of exchange couplings is equivalent to the 
model shown in Fig. W[c). 

Half-doped manganites possess either A, CE or C^Ei-a, type AF magnetic states below 
their Neel temperatures, depending on the identities of the ions A and B in Ai-^B^MnOs 
jlfl| . CE-type AF order consists of zig-zag FM chains antiferromagnetically coupled to 
neig hboring chains. C x Ei_ x order is an incommensurate charge and orbital CE-type order 
llj . CE and C x Ei_ x order is found for wider gap manganites such as Pr .6oCao.4oMn03 
4( and La .5Ca sMn0 3 j3j, while A-type order is found for metallic manganites such as 
Pro.sCao.sMnOa |10J]. Incommensurability and a fine energetic balance between A and CE- 
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type magnetic order are therefore features of the half-doped manganites which also appear 
in ab initio calculations jjj. Magnetic susceptibility data for Pr .6oCao.4oMn03 4] has been 
interpreted in terms of ordering of magnetic moments of Zener polarons in a CE-type state 
below TV = 115 K, and ordering of magnetic moments into Zener polarons at the CO 
temperature, Tco — 330 K. 

The ground state of the classical, magnetic Hamiltonian corresponding to Fig. [IJc) * s 
A-type when the magnitude of Jiis less than half that of J2. The ground state changes to an 
incommensurate spin spiral state with the in-plane component of the wave vector parallel 
to the diagonal direction of J2 when the ratio of exchange constants J2/ |Ji| exceeds 0.5. 
The spin spiral state becomes the commensurate, orthogonal phase described by Efremov 
et al. p when the spiral wavevector becomes (vr/2, tt/2, 7r); the classical magnetic model 
predicts that this state is found only in the limit J2/I Ji | ~~ * 00. Both spin wave theory 
and cluster diagonalization support a picture in which the magnetic structure in half-doped 
manganites is strongly dependent on the relative magnitudes of exchange constants J\ and 
Ji- Specifically, AF coupling between polarons tunes the magnetic correlations between 
a state with in-plane ferromagnetism, similar to the A-type structure, and a state where 
AF correlations are dominant. The latter state is the spin spiral, which is reminiscent of 
the CE-type structure. A similar picture has emerged in Q| with non-bipartite magnetic 
structures competing for the ground state. These authors find that the magnetic structure 
of the ground state around half-filling is subtly dependent on the extent of doping, x, and 
for a part of the phase diagram the ground state is intermediate between a conventional CO 
state and a Zener polaron state. 

A two-dimensional Heisenberg model, with frustrated, AF interactions similar to those in 
the three-dimensional model used here, describes the magnetic properties of organic molec- 



ular crystals 
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131 ] . A Heisenberg model with the same connectivity as that used here has 
been applied to „ -NaV 2 0« Q. A sMla^pin wave approach to the one n S ed in this pape r 
was recently applied to Lao.sCao.sMnOa |15J. In that work each Mn ion was assigned a spin 
3/2 or 2, there were FM couplings within zig-zag chains and AF couplings between spins in 
adjacent chains. Hence that work differs from the present work in the magnetic units used 
(single ions versus dimers) and the presence or absence of frustration in interchain coupling. 

The plan of the rest of this paper is as follows: in section |H] the model is introduced and 
solved at the classical level; in section IIHI linear spin wave theory is applied; in section IIVI 
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the results of diagonalization of the Hamiltonian for small clusters with periodic boundary 
conditions are presented. Finally conclusions are presented in section IVl 



II. HAMILTONIAN AND CLASSICAL SOLUTION 

The Heisenberg Hamiltonian used in this work is 

H = / J J\ ■ S{ + ij t k + ^i,j,k ' Si,j+l,k) + J2 ^i,j,k ' + <^3 s i,j,k ' Sjj./c+l ] . (1) 

The FM zig-zag chains are defined by consecutive steps in the x and — y directions in Fig. 
Wi c )- z is the interplanar axis. i,j and k label sites along the three axes, J\ is FM and 
negative while J2 and J3 are AF and positive. Interactions within xy planes are frustrated 
since J\ prefers parallel spin alignment while J2 favors antiparallel alignment along the 
diagonals. The term bond is used to describe these interactions from here on. Ratios of 
the magnitudes of ab initio exchange constants are J 2 /| J%\ = 0.38 and J 3 /| J\\ = 0.62 (using 
an average value of the FM1 and FM2 of —13 meV for Ji) |5j. The interplanar exchange 
coupling J 3 is not frustrated and the ratio J 3 /\Ji\ will be fixed at 0.5 from now on, except 
where noted. 

At the classical level the angle between neighboring spins is found by minimizing 



the structure factor 



J(0) = Ji (cosq x + cosq y ) + J 2 cos(q x + q y ) + J 3 cosq z . (2) 

The solution is q = (0, 0, n) for J2/I Ji| < 0.5 and q = (q, q, n) with q =arccos(— Ji/2 J 2 ) when 
^2/ |^i I > 0.5. The former case corresponds to FM intraplanar order and AF interplanar 
order; the latter corresponds to a spin spiral in the plane, with the intraplanar component of 
the wavevector along the direction of the J 2 bond. The angle between neighboring spins is q. 
In the limit t/2/l^il ~~ * 00, q — > 7r/2 so that neighboring spins are at right angles; this is the 
orthogonal phase described in [(] . Thus the solution to Eq. ^for small J 2 /| J\\ corresponds 
to A-type AF order, while for J 2 /| J\\ — > 00 the structure is a noncolinear CE-type magnetic 
phase. 

Clusters of cubic symmetry with 2x2x2 and 3x3x3 sites will be considered in 
exact diagonalization calculations in section IIV1 For the second cluster, periodic boundary 
conditions are frustrated for the AF bonds, due to the odd number of spins along the 



corresponding directions. For both clusters the classical ground state is FM in the xy plane 
for J2/\J\\ < 1, with energies per bond J\ and Ji- For J2/I Ji| > 1, the bond energies are 
and — J2 for the 8-site cluster, and Ji/4 and — J2/2 for the 27-site cluster. For the 8-site 
cluster this is orthogonal type of magnetic order [6j, while for the 27-site cluster it is a 
frustrated configuration along the diagonals. Along the z direction the energy per bond is 
— J3 for 8 sites and — J3/2 for 27, regardless of the value of J%I\J\\. Frustration in boundary 
conditions increases the energy for the 27-site cluster. The total energies are 4(2 Ji + J2 — J3) 
for 8 sites and 27(2J X + J 2 — (J 3 /2)) for 27 sites when J2/I Ji| < 1, and the corresponding 
energies for J2/I Ji| > 1 are — 4(J 2 + J3) and (27/2) (Ji — J2 — J3) respectively. The ratio 
J 2 /| Ji| where the ground state correlations change character is different compared with the 
infinite lattice value due to the finite size of the clusters. 



III. LINEAR SPIN WAVE THEORY 



Owing to the spiral nature of the classical ground state for J2/IJ1I > 0.5, it is conve- 
nient to redefine the Sj operators so that the local quantization axis points along the the 
classical solution spin directions. Only one type of bosonic operator is necessary here, al- 
though the Hamiltonian becomes more complex. After introducing the Holstein-Primakoff 
transformation with operators a and Fourier transforming, the Hamiltonian becomes 

H = NE d + SYX A A a k + Y (a H* + a S°-^' (3) 
a- 

where N is the number of sites, E c i = S 2 £ JiCos9i, = — Jj[2 cos9i — (l + cos9i) cos(k- 

i i 

5i)] and = — )] Jj(l — cosOi) cos{k ■ 5i). i refers to bonds in the unit cell with 6i the 

i 

angle at the classical level between two spins connected by bond J i; and 5, the unit vector 
in the bond's direction. After a Bogoliubov transformation to new operators 
diagonalized Hamiltonian is 
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the 



H = NE d + S ( £ cgaUj + £ e 3-*£ ), (4) 



with = JAl — Bi. The average spin magnitude per site along the classical solution 
direction is 



The ground state energy is found by setting the occupation number ce^a^ equal to zero 
for every k. The ground state energy and magnetization per site are plotted in Fig. |2] 
as a function of J2/IJ1I with J 3 /\Ji\ = 0.5. The correction to the energy for quantum 
fluctuations for Sj=7/2 is very small owing to the large spin magnitude; the classical spin 
structure survives with only minor changes in the magnetization. The point of maximal 
frustration where the energy has a maximum is shifted from 0.71 for the classical case to 
0.69 in linear spin wave theory. Quantum fluctuations lower the local magnetization from 
7/2 per site and there is a minimum in the local magnetic moment when the classical ground 
state changes from a FM to a spiral state. 

Spin wave dispersion relations along (k,k,ir) (parallel to the AF J 2 bond), (k,—k,7i) 
(perpendicular to the J2 bond), and (0, 0, k) (perpendicular to planes containing zig-zag 
chains) are shown in Fig. EJ When J2/IJ1I = 0, dispersion relations are characteristic of 
a FM ground state; at the transition point, J2/\J\\ = 0.5, the spin wave velocity vanishes 
along the (k, k, n) direction and magnetization corrections have a local maximum |2l|. When 
J2/\Ji\ > 0.5, changes in the magnetic ground state to a spin spiral become evident in the 
dispersion relation along (k,k,7r) (Fig. Efa)); zero modes appear at the spiraling wave 
vector, whose magnitude increases with Jz/\Ji\- The zero mode occurs at (ir/2, ir/2, 7r) 
when J 2 /| J\\ — > 00. When t/2/l^il < 0.5, spins are ferro magnetically correlated along the x 
and y directions and a dispersion relation characteristic of FM order in the ground state is 
found (Fig. Efb)). But when Jij\J\\ > 0.5, a gap develops at (0,0,0) since there is a spiral in 
the planes. AF order from plane to plane along the (0, 0, z) direction (Fig. Etc)) is reduced 
due to quantum fluctuations as ^/l-^il increases, with the energy gain increasing along the 
J2 bond. Inclusion of the interplanar interaction makes the system three dimensional and 
the integrals in the calculation of the ground state energy and the magnetization per site 



are well-behaved, while in the calculation injl2j and 
enhances the role of quantum fluctuations j22j. 



21( the planar character of the model 



IV. DIAGONALIZATION OF FINITE CLUSTERS 



Hamiltonians (Eq. for three dimensional clusters with periodic boundary conditions 
and spin 7/2 or 1/2 were diaeonalized after block factorization using permutation and spin 

nhn 

symmetries |23L 1241 I25| . Hamiltonians for 2x2x2 and 3x3x3 clusters, which have 



7 



all the symmetries of the bulk, and a 3 x 3 x 2 cluster were diagonalized. Energy eigen- 
value spectra and nearest and next nearest neighbor correlation functions were computed. 
Whenever full diagonalization was not possible, lowest eigenvalues were calculated using 
the numerical package ARPACK 13- The aims were to corroborate results presented 
above by investigating the influence of quantum fluctuations on the classical results, further 
study the competition of FM and AF intraplanar and AF interplanar exchange couplings 
and examine finite-size effects. Diagonalization of Hamiltonians with Sj = 7/2 was limited 
to 8-site clusters due to computer memory requirements, whereas Hamiltonians for 27-site 
clusters could be diagonalized for Sj = 1/2. 

The permutation symmetry group of the Hamiltonian, which is the symmetry group of the 



clusters in real space, is the product of the translation and the point symmetry groups 



24, 



251 ] . The point symmetry group for the cubic clusters is D2h- The Hamiltonian Eq. is also 



symmetric with respect to time reversal. The spin inversion symmetry group consists of the 
identity and the spin inversion operation, and its product with the permutation group gives 
the total symmetry group of the lattice. Eigenstates are characterized by the momentum k 
and the irreducible representation of the total symmetry group that they belong to. 

The ground state energy for the 2x2x2 s$ = 7/2 cluster is plotted in Fig. Ufa) 
for two values of J^/\J\\ as a function of ^/l-Al- It is a total spin S = state with 
momentum k = and is symmetric with respect to spin inversion. Frustration is maximal 
for J2/IJ1I = 0.91592 and 0.91895 when Jz/\Ji\ = 0.4 and 0.5, respectively. The classical 
transition point J2/IJ1I = 1 is therefore renormalized by quantum fluctuations towards 
smaller J2/\J\\- 

Around the points of maximal frustration the characteristics of the ground state change, 
as shown by the correlation functions in Fig. EJb). For small J2/I Ji| the nearest neighbor 
and next nearest neighbor correlation functions are strongly FM. For larger J2/I Ji| the near- 
est neighbor correlation function is very small, while the next nearest neighbor correlation 
function becomes strongly AF. At the same time, the interplanar correlation function, which 
is always AF, loses some of its strength due to quantum fluctuations, with the spins gaining 
energy predominantly via the AF J 2 bond. The changes in the nature of the correlation 
functions are quite sharp, similar to the behavior of the energy around maximum frustration. 
The relative changes of the correlation functions as a function of ^/l^il are very similar to 
changes in the classical solution. This classical-like behavior is due to the large value of the 
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spin, Si = 7/2. 

Calculation of the low energy properties for the 27-site, Sj = 7/2 cluster is not possible, 
as noted above. To study the effect of the frustrating bond J2 as a function of cluster size, 
we consider s, = 1/2 and estimate the changes compared to the 8-site case. The ground 
state energy for J%/\J\\ = 0.4 and 0.5 is shown in Fig. Efa) for 8 sites. It is an S = 
state with momentum k = and is symmetric with respect to spin inversion. The points of 
maximal frustration are now around 0.57 and 0.58 respectively, significantly changed from 
the classical and = 7/2 values. The correlation functions are plotted in Fig. [31(b), and the 
change of the in-plane correlation functions from a FM to a diagonal AF character is now 
smoother. This is attributed to the smaller magnitude of the spins. The gain in energy via 
the J 2 bond is now more significant. 

For the 27-site cluster the ground state energy for J%/\J\\ = 0.4 and 0.5 is plotted in 
Fig. infa) as a function of Ja/I^il- Frustration is maximal for t/2/l^il — 0.7212 and 0.7278. 
The ground state energy has momentum k = (0, 0, 7r), which is doubly degenerate. Its point 
symmetry grou p is C^h, a subgroup of D 2 h- The ground state belongs to the A g irreducible 
representation [2j|. Correlation functions for the ground state are plotted in FiglHfb). Their 
behavior is similar to the correlations in the 8-site cluster and they change character around 
the point of maximal frustration. The nearest neighbor intraplanar correlation function 
is 0.247 and almost fully polarized at ^/l^il = 0, and then drops to approximately one 
third of this value for larger J2/I The diagonal correlation function starts with the same 
strong FM character at J2/\Ji\ = and reverses sign for larger values. It is equal to —0.205 
for J2/IJ1I = 1, which shows strong AF correlation when compared with the value for 
Ji = 0, which is —0.238. The interplanar correlation function originally decreases slightly 
with increasing t/2/l^il and for higher values increases slightly to accommodate the increase 
of the diagonal intraplanar correlation function. Closer inspection of the plot (Fig. |7J) 
reveals two discontinuities in the correlation functions for values of J2/\Ji\ equal to 0.736341 
and 0.75008, where the two lowest energy states change roles as the ground and the first 
excited state. The two states have the same momentum and belong to the same irreducible 
representation, but the ground state between the two discontinuities is an S = 3/2 state, 
while the other is S = 1/2. At the same points the derivative of the energy with respect to 
J2/I Ji\ is discontinuous. 

As was the case for the 8-site cluster, the transition from a FM to a spin spiral state has 
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been renormalized by quantum fluctuations, however the corresponding J2/I Ji| value is closer 
to the classical value of 1 compared with the 8-site cluster. The changes of the intraplanar 
correlation functions as functions of ^/l^il are sharper compared with the changes in the 
8-site cluster. It is expected that the results for the Si = 7/2 case will be similar and the 
point where correlation functions change will be closer to 1, compared with the 8-site case. 
This result agrees with linear spin wave theory, where quantum fluctuations have a small 
effect on the classical solution. 

When Si = 1/2 full diagonalization is possible for systems of 18 sites; 3x3x2 cluster 
Hamiltonians were diagonalized. The specific heat of the cluster is plotted in Fig. IBf a) as a 
function of temperature for several values of J2/I J\\. There is a shoulder at low energy which 
disappears for ^/l^il > 0.7, again indicating a change in the nature of the ground state. At 
the same time, the main peak is pushed towards lower temperature and its value decreases. 
The magnetic susceptibility as a function of temperature is shown in Fig. IHtb). The shoulder 
found for lower J2/\J\\ values disappears for higher values and the peak position shifts to 
lower temperature, signifying again a qualitative change as the frustrating interaction gets 
stronger. 



V. CONCLUSIONS 



A Heisenberg model with FM nearest neighbor interactions J\ and AF next nearest 
and interplanar interactions J2 and J3 has been studied as a prototype for the magnetic 
behavior of polarons forming in the half-doped lanthanum manganite Lao.sCao.sMnOs j^. 
The polarons are spin 7/2 objects formed by two Mn and an O ion. At the classical level the 
ground state is of the A-type for t/2/l^il — 0.5, while for higher J 2 it is a spin spiral. When 
I > 0.5, there is a spiral in the xy plane with every other pair of spins parallel in zig- 
zag chains, and an angle q = arccos(—Ji/2J 2 ) between neighboring spins. For J2/I Ji| — > 00 
the spiral becomes the orthogonal state described by Efremov et al. in |(|, where polaron 
moments are perpendicular to each other along the zig-zag chains and antiparallel along 
the direction of j \. The effect of quantum fluctuations on spin wave theory is small due 
to the large magnitude of the spins, Si = 7/2. Hence the magnetization per site is not 
significantly changed from its classical value. Diagonalization of finite clusters also shows 
that quantum fluctuations do not significantly alter the classical solution, and the ratio 
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Jij\ J\\ for which there is a transition from FM order in the planes to one where spins are 
coupled antiferromagnetically via the J2 bond is close to its value at the classical level. 

The calculations in this paper show that the magnetic structure of Zener polarons in 
Lao.sCao.sMnOa is subtly dependent on the ratio ^/l^il- The value from ab initio calcu- 
lations was found to be 0.38 15J], however UHF calculations underestimate the value of AF 



couplings, and recent calculations of real space structure tend to favor a ratio close to 0.5 
with a CE-type magnetic ground state [7]. Thus the physically relevant parameter space of 
the Hamiltonian (Eq. |TJ) has J2/IJ1I ~ 0.5, where for Jz/\Ji\ > 0.5 the magnetic structure 
is non-bipartite and depends on the exact value of J2/IJ1I. Similar conclusions have been 
drawn by Efremov et al. in 0. In that paper the authors have found that around half-doping 
magnetic order depends sensitively on the extent of doping x, and the orthogonal state is 
the ground state for a part of the phase diagram, while a state which is a superposition of 
the conventional CE-type order and the orthogonal phase is lowest in energy for different 
combinations of parameters. The angle between neighboring spins is 2ir/3 in that phase. 
In our model the angle between spins linked via the J2 bond and belonging to different 
zig-zag chains is 2tt/3 when J2/I J\\ = 1, while the angle between nearest neighbors is 7r/3 in 
that case. This phase is also weakly renormalized by quantum fluctuations. Therefore the 
Hamiltonian (Eq. ^) predicts phases which differ from the conventional CE-type order. 

Diagonalization of small clusters with cubic symmetry also showed that there is a change 
in the character of the ground state as a function of J2/\J\\- For small ratios spins are fer- 
romagnetically ordered within planes, while for higher values the next-nearest neighbor AF 
interaction dominates. As was the case with linear spin wave theory, the diagonalizations 
show that the classical results are not significantly altered by quantum fluctuations. A spin 
magnitude of s; = 1/2 was also considered to study finite size effects, and the transition from 
the 8- to the 27-site cluster showed that the ratio where the character of the ground state 
changes comes closer to the classical value. Results from larger clusters would be needed to 
firmly establish this point, however memory requirements prohibit diagonalizations of larger 
cubic clusters. The role of the AF interaction J2 was also evident in full diagonalizations 
of the model for systems of 18 sites with Sj = 1/2. Specific heat and magnetic suscepti- 
bility data were calculated, and the graphs change qualitatively as a function of ^/l^il) 
showing again the role of the diagonal AF interaction in the development of AF intraplane 
correlations. 
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The spin wave dispersion was calculated for the material Ndo.4sSro.55Mn03 in [29], and it 
was fitted with a Heisenberg model with FM nearest neighbor interactions in the planes and 
AF interactions between planes. There was also an anisotropy term. The relative strength 
of the AF with respect to the FM interaction was found to be 0.620, which is in agreement 
with the values in [5j]. The anisotropic interaction was relatively small. 

The authors wish to acknowledge discussions with R. Kenna. Calculations were carried 
out at the Trinity Centre for High Performance Computing. NPK was supported by a Marie 
Curie Fellowship of the European Community program Development Host Fellowship under 
contract number HPMD-CT-2000-00048. 
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FIG. 1: Schematic illustration of the connectivity of the model: (a) effective exchange interactions 
between polarons (FM stands for ferromagnetic and AF for antiferromagnetic) . (b) Magnetic unit 
cell, (c) Connectivity of the model in the planes. Interplanar interactions are AF. 
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FIG. 2: Ground state energy and magnetization per site as a function of J2/\J\ \ for Jz/\J\\ = 0.5 
within linear spin wave theory. Top: solid line: classical energy, dotted line: spin wave energy. 
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FIG. 3: Dispersion relations along different directions in the Brillouin zone for J%/\J\\ = 0.5 
within linear spin wave theory. Top: (k, k, it), middle: (k, —k, tt), bottom: (0,0, k). straight line: 
I = for top and < Jil\J\\ < 0.5 for middle and bottom, dotted line: J2/IJ1I = 0.5 for 
top, dashed line: J2/\Ji\ = 1, long-dashed line: J2/IJ1I = 5. 
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FIG. 4: Energy and nearest and next nearest neighbor correlation functions of the ground state 
for Si = \ and a 2 x 2 x 2 cluster as a function of J2/IJ1I. Top: J%/\Ji\ = 0.4 (o) and J3/IJ1I = 0.5 
(□). Bottom: correlation function for bond J\ (o), for bond J2 (□), and for bond J3 (o). The lines 
are guides for the eye. 
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FIG. 5: Energy and nearest and next nearest neighbor correlation functions of the ground state 
for Si = \ and a 2 x 2 x 2 cluster as a function of Jil\J\\- Top: J?,/\Ji \ = 0.4 (o) and J3/IJ1I = 0.5 
(□). Bottom: correlation function for bond J\ (o), for bond J2 (□), and for bond J3 (o). The lines 
are guides for the eye. 
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FIG. 6: Energy and nearest and next nearest neighbor correlation functions of the ground state 
for Si = \ and a 3 x 3 x 3 cluster as a function of Jil\J\\- Top: Jz/\J\ \ = 0.4 (o) and Jz/\J\ \ = 0.5 
(□). Bottom: correlation function for bond J\ (o), for bond J2 (□), and for bond J3 (o). The lines 
are guides for the eye. 
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FIG. 7: Nearest and next nearest neighbor correlation functions of the ground state for S{ = \ and 
a 3 x 3 x 3 cluster as a function of 3ij\J\\ for Jzj\J\\ = 0.5. o : correlation function for bond Ji, 
□ : for bond J2, o : for bond J3. The lines are guides for the eye. 
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FIG. 8: Specific heat (top) and susceptibility (bottom) as function of temperature for J3/I J\\ = 0.5. 
solid line : J2/\Ji\ = 0.4, dotted line : J2/\Ji\ = 0.5, dashed line : J2/IJ1I = 0.6, long-dashed line 
: J 2 /\Ji\ = 0.7, dot-dashed line : J 2 /\Ji\ = 0.75, dot-long dashed line : J 2 /\Ji\ = 0.8. 
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